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Entanglement polytopes result in finitely many types of entanglement that can be detected by only measuring 
single-particle spectra. With high probability, however, the local spectra lie in more than one polytope, hence 
providing no information about the entanglement type. To overcome this problem, we propose to additionally 
use local filters. We experimentally demonstrate the detection of entanglement polytopes in a four-qubit system. 

Using local filters we can distinguish the entanglement type of states with the same single particle spectra, but 
which belong to different polytopes. 

PACS numbers: 03.65.Ud, 03.67.Mn, 71.10.Pm, 73.43.Nq 


Entanglement among quantum systems is a kind of quan¬ 
tum correlation that is stronger than any possible classical 
correlation 111]. At the fundamental level, entanglement is the 
very mystery of quantum mechanics; and at the practical level, 
it can be used as a physical resource to perform computation 
and communication tasks that are impossible for classical sys¬ 
tems. Essentially, entanglement comes from the tensor prod¬ 
uct structure of the Hilbert space of N systems-qubits in our 
case. An A-qubit quantum state I'Eat) is entangled if it can 
not be factored into products of quantum states of each of the 
qubits. 

One central question regarding entanglement is that how 
I'Eat) may be entangled and how to detect that feature in prac¬ 
tice. An obvious fact is that the parameters needed to specify 
I'Eat) grows exponentially with N. A natural idea to eliminate 
some of the free parameters is to take that two states I'Eat) 
and |<I>jv) have similar entanglement features if they can be 
connected by some single-qubit operations, for instance, local 
unitary (LU) transformation or stochastic local opera¬ 
tion combined with classical communication (SLOCC) IIsI-ItI]. 

Eor A = 2, the Schmidt decomposition tells us that |'E 2 ) = 
A<)|00) + All 11) up to LU, with Aq > Ai and Aq+A i^ = 1. 
Different Aq € [l/v/2,1] corresponds to different LU classes 
of entanglement, which are in fact infinitely many. Up to 
SLOCC, however, there is only one class of entangled states 
which contains the EPR pair with Aq = 1 /v/2- Lor A = 3, 
up to SLOCC, there are only two types of entanglement: the 
GHZ-type state and the W-type state fit]. These two types 
can be distinguished by a quantity called 3-tangle, which is 
however not a single-copy observable and hence cannot be 
directly measured in experiment fl] (that is, to get the value 
of 3-tangle, one either needs to measure jointly on multiple 
copies of the states, or needs a state tomography j^). 

Lor any A > 3, SLOCC no longer results in a finite number 


of entanglement types fin. Despite the efforts of study¬ 
ing SLOCC classification of entanglement for A > 3 systems, 
the exponential growth of parameters with A for describing 
I'La^) makes it hopeless to extract clear physical meanings of 
these classifications. It is highly desired to coarse-grain these 
classes such that we can grasp the key features of each en¬ 
tanglement type. The concept of entanglement polytopes pro¬ 
vides an elegant idea to meet this need where for each 
A there exists only finite number of types. More importantly, 
the polytopes are directly detectable in experiments via mea¬ 
suring only single-particle spectra of each qubit UMl- 

In this work, we experimentally demonstrate the detection 
of entanglement polytopes in a four-qubit system. Unfortu¬ 
nately, it turns out that different entanglement polytopes form 
a nested hierarchy llol - [l^ . and they may have a large overlap. 
If the vector of local spectra of a state I'Tat) lies in an overlap¬ 
ping region, then we cannot uniquely identify the polytope 
that I'Tiv) belongs to (see e.g. a recent experiment in which 
the states are chosen to be in non-overlapping regions S). 


It turns out that for a randomly chosen three-qubit pure state 
|'T 3 ), the probability that the vector of local spectra of |'T 3 ) 
lies in the overlapping region of the W and GHZ polytopes is 
« 94%. In general, for a randomly chosen state A-qubit state 
I'Ta:), with high probability the vector of local spectra falls 
in some overlapping region of polytopes (we include a more 
detailed discussion of these probabilities in the Appendix). To 
overcome this difficulty, we use local filters (see, e.g. lll4l - ll8n i 
to effectively distinguish states with the same single particle 
spectra, but which belong to different polytopes. 


Entanglement Polytopes - An A-qubit quantum state |'Ta:) 
is said to be convertible to another A-qubit state |<I>Ar) via 
SLOCC if there exists a sequence of local operations and 
classical communication that converts the state I'Lat) to |<I>a:) 
with nonzero probability. The states I'La^) and |<I>a^) are said 
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to be SLOCC equivalent if |<I>aj) is convertible to |<I>Ar) and 
vice versa. It has been shown that I'Pat) is SLOCC equivalent 
to |<I>Ar) if and only if there exist invertible matrices M, for 
i= such that 


= Ml 0 M 2 0 ... 0A/Ar|<I>/jr). (1) 

This SLOCC equivalence relation partitions all A^-qubit pure 
states into SLOCC equivalent classes, called the SLOCC or¬ 
bits. 

For an A^-qubit state |'Faj), each of the single-particle re¬ 
duced density matrices p,- for i = has two eigenval¬ 
ues that are normalized, i.e. A“+Af = 1. It suf¬ 

fices to consider the maximum eigenvalue of p,-, i.e. = 
max(A“,A,^), with ^ < A™ < 1. The A^-dimensional vector 
A = (A™, A™“,..., A™®’‘) then corresponds to a point in 

It has been shown that for the closure of an A^-qubit SLOCC 
orbit 0 of'Fiv (i e. someM, in ([T]i possibly non-invertible), all 
the points A of all (8>iLiM, |‘FAr) € 0 form a polytope in R^, 
called the entanglement polytope of 0. Moreover, for any fi¬ 
nite N, there are only finitely many polytopes. This provides 
a natural classification of entanglement for A^-qubit states, 
which ‘coarse-grains’ the infinitely many SLOCC orbits (for 
N > 3). Since the entanglement polytope for an SLOCC orbit 
is fully determined by the local spectra of the states in the or¬ 
bit, this offers an appealing experimental approach for identi¬ 
fying the entanglement type for an A^-qubit system, for which 
only measurements on single particles are required. 

As a simple example, for N = 3, with only two kinds 
of SLOCC orbits (GHZ and VT-type states), there are two 
polytopes in R^. The polytope corresponds to the 

GHZ-type states, with vertices (1 /2,1 /2,1 /2), (1 /2,1 /2,1), 
(1/2,1,1/2), (1,1/2,1/2), (1,1,1); and corresponds 
to the W-type states, with vertices (1,1,1), (1/2,1/2,1), 
(1/2,1,1 /2), (1,1 /2,1 /2). Obviously, C 
shows that entanglement polytopes for different SLOCC or¬ 
bits may overlap. 

In practice, for a state |'Fa'), while a point A may clearly 
distinguish its entanglement type, if A is in an overlapping re¬ 
gion of two polytopes, we fail to get information on which 
entanglement type the state belongs to. In the N = 3 case, 
for instance, this means that a point A G 0^^ fails to distin¬ 
guish the fT-type entanglement from the GHZ-type entangle¬ 
ment. Unfortunately, for a randomly chosen pure state of three 
qubits, with Ri 94% probability the corresponding A falls into 

Luckily, one can apply local filter operations M, to the 
system to ‘move around’ A, with the hope that A ends up in 
a non-overlapping area of polytopes. As demonstrated in the 
three-qubit case, this step becomes crucial in practice when 
using the polytope method for detecting entanglement types, 
as the probability of overlapping is high. 

Four-qubit Polytopes - Our experiments demonstrate the 
detection of entanglement types of four-qubit states. In this 
case, there are infinitely many SLOCC orbits. The full poly¬ 


tope, containing A for any four-qubit state, denoted by 
is spanned by the vertices 


( 1 , 1 , 1 , 1 ),( 1 , 1 , 1 , 1 ),( 1 , 1 , 1 , 1 ), 


( 1,1 1 , 1 ),( 1 , 1,1 1),(1 1 , 1 , 1 ),( 1 , 1 , 1 , 1 ). 


Up to permutation of the qubits, there are 6 other polytopes 
inside which may also be mutually overlapping. Sim¬ 

ilar as in the N = 3 case, for a randomly chosen four-qubit 
state, the chance that A lies in an overlapping region is high 
(for details, see the Appendix). Therefore, we have to apply 
local operations to ‘move around’ A. The proposed experi¬ 
ment is given by the diagram in Fig.[T] 



1 (post-selection) 


FIG. 1: Circuit diagram of the experimental setup. 


Here tJ,- for i = 1,2 denotes a unitary local transformation 
U of the form 

_ /cost?,' — sint?,'\ /l 0 \ / cost?,' sint?,'\ 
l^sint?,' cost?,' y 1^0 — ly \^—sint?,' cost?,y ’ 

and 7 denotes a non-unitary local transformation 



In Fig. [T] two of the qubits encounter non-unitary local 
transformations: qubit 1 is measured in some basis and post- 
selected, resulting in A™” = 1; qubit 4 is going through a filter 
operation given by A^. In the most general case, one can also 
apply local filter operations (or measurements) on the other 
qubits. However, a single filter (or measurement) may already 
suffice to ‘move around’ A to non-overlapping regions of the 
poly topes, depending on the input state 

Experimental setup - In our experiments, two different 
four-qubit states and are prepared, where 

|vj/{i)^ = ::11{\hhhh) + \vvvv)) 

+ ^{\HV)0\yHmHV) + \VH)), (2) 

o 
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FIG. 2: Detailed configurations for preparing the states and and for realizing the operator Ay are shown in the orange boxes and the 
blue box, where HWP3, HWP4 and HWP5 are rotated by ^(arcsiny), j, and ^ respectively. The unitary operator is realized by HWP2 
at specific angles. A QWP and a HWP in front of a polarization beam displacer (PBS) in each mode are used to implement the measurement 
in different bases for the standard state tomography. Post-selection in some basis of qubit 1 (the state of the photon in mode ‘3’) is realized by 
collecting only photons in one of the output modes of the PBS. The indices in the figure denote the spatial modes. 


and I'PP)) is the four-qubit GHZ state 
\/2 

|'p(2)) = + (3) 

The qubits are encoded by horizontal \H) and vertical |T) po¬ 
larization. The goal is to determine the entanglement type for 
each of the input state using the polytope method. For both 
and we have A = That is, local 

spectra do not tell them apart, hence local filter operations are 
needed to ‘move around’ A. 

The SLOCC orbit of the four-qubit GHZ state cor¬ 
respond to the full polytope However, the state 

corresponds to a smaller polytope C with vertices 

(ll,l,l),(ll,il),(i 1,1,1). 

The smaller polytope is characterized by the additional 
constraint 

/(A) = -Ar“ -F A2”“ + Ag'"” + Af “ > 1 (4) 

and all permutations of it. 

Our experimental setup for the states and |'FP)) 

is shown in Fig. |2l A 390 nm femto-second pump light. 


frequency-doubled from a 780 nm mode-locked Ti:sapphire 
pulsed laser (with the pulse width about 150 fs and repeti¬ 
tion rate 76 MHz) was used to pump the respective down- 
converter. For the preparation of | ^ ^, a 2 mm type-II phase- 

matched BBO crystal is used as down-converter to produce 
two pairs of entangled photons ifl^ . and two 1 mm BBO 
crystals are used to compensate the birefringence of o-light 
and e-light in the 2 mm BBO. HWPl rotates the polarization 
of the photons in path ‘2’ (horizontal to vertical and vertical 
to horizontal). Then after the beam splitters (BS), the above 
two pairs of entangled photons are transformed into the state 
In mode ‘6’, we use two beam displacers and three 
half wave plates (HWPs, HWP3 is used for balancing the op¬ 
tical length of the two beams between partdisplacers) to con¬ 
struct the local filter Ay. For the four-qubit GHZ state 
shown in the right part of F ig. El a cascaded sandwich beam¬ 
like BBO entangled source 1201] is used. A PBS combines the 
photons from mode ‘1’ and ‘2’. We will get the four-qubit 
GHZ state jf there is one photon in each of the modes 

‘3’, ‘4’, ‘5’, and ‘6’ il. 

Results - We first perform full quantum state tomography to 
reconstruct the density matrix of and |'P^2))^ the fidelity 
of which are 0.9422 ± 0.0036 and 0.9001 ± 0.0038. Then 
we collect data from each mode to obtain the corresponding 
single-qubit density matrix and calculate their local spectra 
for both states. As shown in Table |I] the local spectra of 
and are almost identical, so we can not distinguish their 

entanglement polytopes. 

To distinguish the entanglement polytopes of |'P4)) and 
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state 

^max 


5 max 

^max 

fW 

|VI/(1)) 

0.529(4) 

0.514(4) 

0.540(4) 

0.530(4) 

1.056(8) 

|vp(2)) 

0.521(4) 

0.524(4) 

0.535(4) 

0.525(4) 

1.062(8) 


TABLE I: The local spectra A™”, Aj"”, A™’' together with 

/(A) for the states and The uncertainties inside the brack¬ 
ets are obtained by Monte Carlo simulation (1000 steps). 




r 



A max 

fW 

a 

Tt/S 

\/y/2 

0.609(10) 

0.831(9) 

0.701(10) 

1.141(18) 

b 

Tt/S 

l/\/3 

0.557(10) 

0.850(9) 

0.614(10) 

1.021(17) 

c 

Tt/S 

l/\/5 

0.603(10) 

0.875(9) 

0.553(10) 

1.032(18) 

d 

3n/32 

1/a/5 

0.713(9) 

0.883(9) 

0.717(8) 

1.313(17) 

e 

0 

1 

0.657(9) 

0.848(8) 

0.857(8) 

1.362(17) 

f 

0 

1 

0.525(9) 

0.544(8) 

0.516(8) 

0.584(19) 


TABLE II: Setting of the parameters di and y for the data points 
labeled ‘a ~ /’, together with the measured local spectra A^“, A™’', 
A™“ and the resulting value of /(A). The uncertainties inside the 
brackets are obtained by Monte Carlo simulation (1000 steps). 


we then try to move A out of the smaller poly tope 
using local filters, as illustrated in Fig.[T] We fix t92 = —7r/8, 
and then measure the first qubit in the computational basis. By 
post-selection we have A[”“ = 1. For each setting of t>i and y, 
we perform tomography of the qubits 2, 3, and 4 to determine 
the values of A™“, A““, and A““ (see Tablelllli. The smaller 
poly tope is characterized by /(A) > 1. 


1 


imax 

>14 


0.5 


e 

• d 

• a 

f-.- 

• f 
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FIG. 3: Three-dimensional polytopes. The pink region and the blue 
region represent the polytope of and respectively. Ex¬ 
perimental data ‘a ~ e' is for |'I'F)) while ‘/’ is for Error 

bars are too small to identify (see Table Hit. 

The results are illustrated in Fig. [3] and |4] In 
Fig. m the data is shown in the three-dimensional poly¬ 
tope for A™’‘, A™®’‘, A™’‘ as by post-selection of the 
first qubit, A™’‘ = 1. The smaller polytope be- 



FIG. 4: Experimental and theoretical results for /(A) as a function 
of 1)1 and 1/y^. (a) shows the plot for the state |*I'F)) together with 
the experimental data points ‘a ~ e’. (b) shows the plot for the state 
|*1'(^)) and the data point ’/’ from the experiment. Error bars are too 
small to identify (see Table Hit. 


comes a three-dimensional polytope with vertices 

(1,1/2,1/2), (1/2,1,1/2), (1/2,1/2,1), (1,1,1), and the full 
polytope becomes a three-dimensional polytope 

wjjjcjj contains and has an additional vertex 
(1 /2,1 /2,1 /2). The data point / of the state outside 

of shows that is not in In contrast, the data 

points a,h,c,rf,e obtained from all lie in which in¬ 
dicates that belongs to This shows that and 

have different entanglement types. 

In FiglU the plot of /(A) as a function of t?i and 1/y^ is 
shown. For the smaller poly tope the inequality /(A) > 1 
always holds. A violation of this inequality signals that the 
state (point /) is not in In contrast, the data points 

a,b,c,d,e obtained for all satisfy /(A) > 1, which in¬ 
dicates that belongs to This shows that and 

have different entanglement types. 

Since the first photon of the four-qubit state is post-selected 
and the last photon goes through a non-unitary filter, the prob¬ 
ability of succes for the experiment is 0.2917,0.2222,0.1667, 
0.1768, 0.5, and 0.5 for our experimental data ‘a ~/,’ respec¬ 
tively. 

Because the birefringence of the o-light and the e-light 











































5 


in the BBO (down-converter) cannot be compensated com¬ 
pletely, and because of the high-term noise from the SPDC 
process and some mode mismatch, we do not obtain the pure 
states and but some noisy version of them. 

Nonetheless, the errors in our experiment are mainly due to 
the time uncertainty of the photon pairs generated in the BBO. 
The coincidence counts obey a Poisson distribution, the pa¬ 
rameters of which we estimate from the experimental data. 
Then we perform Monte Carlo simulation to estimate the er¬ 
rors indicated in Tables HI and HU 

Summary - We experimentally demonstrate the detection of 
entanglement polytopes in a four-qubit system. We use local 
filters to effectively distinguish states with the same single¬ 
particle spectra, but which belong to different poly topes. This 
provides a new tool to experimental detection of entanglement 
in a multi-qubit system using only local operations. 
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APPENDIX 


A. Entanglement polytopes 

Two quantum states | xj/i ) and 1 1 // 2 ) are said to be equivalent with respect to SLOCC if there exists a sequence of local operations 
and classical communication that converts the state | \j/i ) into 11 // 2 ) with non-zero probability p 1^2 > 0 , and another protocol for 
the conversion of \y/ 2 ) into that succeeds with probability p 2 ^-i > 0. As we only require the success probabilities to be 
non-zero, it is sufficient to consider one branch of the protocol that has non-zero success probability. Thus we can, for example, 
write 




(5) 


where the matrices M, correspond to the combination of all operations performed on particle i. In 11221] it was shown that 
the matrices Mi in Q can be replaced by invertible matrices. Moreover, all matrices can be chosen to have determinant one, 
combining all scalar factors with the success probability. Hence we have that the states |v/i) and I 1 // 2 ) are in the same SLOCC 
class if and only if there is a non-zero constant A G C and matrices A, G SL(^f,) (where dt denotes the dimension of subsystem /) 
such that 


A|v/ 2) =Ai(g)A2(g)---(8)A„|v/i). 


(6) 


Results on the classification of pure four-qubit states with respect to SLOCC can be found in Il23l - l26l] . Note that in the literature, 
sometimes the scaling factor A is incorrectly ignored. 

Unlikely the situation for local unitary transformations, polynomial invariants of the group SL(c/)®" yield only a necessary 
condition for SLOCC equivalence of two quantum states. 


Proposition 1 Let /i, ... ,/m be homogeneous polynomial invariants of the group SL{d)^". If the normalized states 1 1 / 1 ) and 
I 1 // 2 ) are in the same SLOCC class, then there exists a non-zero constant A G C such that 

= fori=\,...,m. (7) 


In the case of four qubits, we have four polynomial invariants Boooo> AJqqqq, Oqqqq, and /pooo of degree 2,4, 4, and 6 , respectively, 
see 1271. 

In order to get a necessary and sufficient criterion to decide SLOCC equivalence, one may consider covariants. Two vectors 
are in the same orbit of the group SL(c/)®" if and only if all covariants agree. Again, one has to take care of the scaling parameter 
A to apply this criterion. In the case of four qubits, there are 170 covariants, see 12811 . 

It has been shown that the points corresponding to the (sorted) spectra of the single-particle reduced density matrices of pure 
quantum states in the closure of an orbit under SLOCC transformations form a so-called entanglement polytope, see i29ll30l] . 
The vertices of the polytope correspond to the covariants that do not vanish. 


B. Four-qubit Polytopes 

In the case of four qubits, there are 7 different 4-dimensional polytopes up to permutation of the qubits iIt^ . Lower-dimensional 
polytopes correspond to states that partially factorize. When we also consider the permutations, the polytopes ,^3 and come 
in 6 different versions, while the polytopes and ^^‘2 split into 4 different subtypes. The vertices of all polytopes are listed in 
Table Uni The poly tope is contained in all other poly topes, and the polytope ,^7 = is the largest polytope containing 
all states. Fig. |5] illustrates how the polytopes are contained in each other. 

In the experiment, we investigate the four-qubit state 

|vf(i)) = + \VVVV)) + ^{\HV) -G \VH)){\HV) + \VH)), ( 8 ) 

3 6 

and the four-qubit GHZ state 


I'pC)) = ::^{]^hhhh) + \vvvv)). 


(9) 


The qubits are encoded by horizontal \H) and vertical |U) polarization. Evaluating the covariants from 12811 for the states 
and we find that the corresponding polytopes are = £^4 and the full polytope respectively. The polytope 

£^4 is obtained from £^2 by removing the vertex (1 /2,1 jl, 1 /2,1) and all its permutations. The discriminating inequalities are 


/(A) = -Ai”” -G Aj™-G Ag"'” -G A 4 ™ > 1 


( 10 ) 
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FIG. 5; Lattice of the different entanglement polytopes for four qubits. Note that, in general, the polytopes intersect non-trivially. 

and all its permutations. 

C. Volume of the Polytopes 

In the case of three qubits, we have only two three-dimensional polytopes C corresponding to the SLOCC class 

containing the VF-state and the G//Z-state, respectively. Picking a pure three-qubit state with respect to the Haar measure at 
random, the resulting distribution of the eigenvalues of the local density matrices has been computed in El. From this one finds 
that the volume of the sub-polytope is 203/216 ~ 93.98%. Hence the probability for a random three-qubit state to have a 
local spectra corresponding to a point outside the polytope is only 13/216 ~ 6.02%. 

For four qubits, we computed the local spectra of 10^ random pure states and determined which of the poly topes contains the 
vector of local spectra. The results are summarized in Table HVl While the poly tope corresponding to the state of 

our experiment is fairly low in the hierarchy of polytopes (see Fig. |5]l, the local spectra of only 9522 out of one million random 
states violate the discriminating inequalities (fTOl i. Hence, the chance for a random four-qubit state to have a local spectrum that 
lies outside of is only about 0.95%. This clearly indicates that one has to apply local filters in order to get information about 
the entanglement polytopes. 

Note that after measuring one of the qubits and post-selection of the measurement outcome, we have a four-qubit state that 
factors into a single qubit and a three-qubit state. The polytope !^a, is mapped to the three-qubit polytope which has a 
volume of about 94%. Hence the local measurement increases the chance for a random state to lie outside the smaller polytope 
from less than one percent to about six percent. 
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.9>‘{ 

(1,1,1,1). (1/2,1/2,1,!), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 

(1/2,1/2,1,1/2), (1/2,1,1/2,1/2), (1,1/2,1/2,1/2), (1/2,1/2,1/2,3/4) 

gib 

(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 
(l/2,l/2,l,l/2), (1/2,1,l/2,l/2), (1,1/2,1/2,1/2), (1/2,1/2,3/4,1/2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 
(l/2,l/2,l,l/2), (1/2,1,1/2,1/2), (1,1/2,1/2,1/2), (1/2,3/4,1/2,1/2) 

g,d 

(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 
(l/2,l/2,l,l/2), (1/2,1,l/2,l/2), (1,1/2,1/2,1/2), (3/4,1/2,1/2,1/2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1,1/2), 

(1/2,1,1/2,1/2), (1,1/2,1/2,1/2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 

(1/2,1,1/2,1/2), (1,1/2,1/2,1/2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 

(1/2,1/2,1,1/2), (1,1/2,1/2,1/2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 
(l/2,l/2,l,l/2), (1/2,1,l/2,l/2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1,1/2,1/2), 

(1,1/2,1/2,1/2) 

^ 3 * 

(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1,1/2), 

(1,1/2,1/2,1/2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1,1/2), 

(1/2,1,1/2,1/2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 

(1,1/2,1/2,1/2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 

(1/2,1,1/2,1/2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 

(1/2,1/2,1,1/2) 

^4 

(1,1,1,1), (l/2,!/2,l,l), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1/2) 

^5 

(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2) 

■n 

(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1,1/2,1/2), 

(1,1/2,1/2,1 /2), (1/2,1 /2,1 /2,1 /2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1,1/2), 

(1,1/2,1/2,1/2), (1/2,1/2,1/2,1/2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1,1/2), 

(1/2,1,1/2,1/2), (1/2,1 /2,1 /2,1 /2) 


(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 

(1,1/2,1/2,1 /2), (1/2,1 /2,1 /2,1 /2) 

^6 

(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 

(1/2,1,1/2,1 /2), (1/2,1 /2,1 /2,1 /2) 

^6 

(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 

(1/2,1/2,1,1 /2), (1/2,1 /2,1 /2,1 /2) 

&>1 

(1,1,1,1), (1/2,1/2,1,1), (1/2,1,1/2,1), (1/2,1,1,1/2), (1,1/2,1/2,1), (1,1/2,1,1/2), (1,1,1/2,1/2), (1/2,1/2,1/2,1), 

(1/2,1/2,1,1/2), (1/2,1,1/2,1/2), (1,1/2,1/2,1/2), (1/2,1/2,1/2,1/2) 


TABLE III: The different entanglement polytopes for four qubits given by their vertices 
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996761 

990140 

.9>\ 

990137 

990204 

990262 

.^2 

863481 

705172 

704928 

704932 

704791 

^3 

781562 

607121 

607010 

607176 

606791 

606925 

&>{ 

607051 

990478 

990478 

130165 

130165 

■^6 

1000000 

n 

995287 

995277 

995320 

995158 

^1 

995201 

995191 

1000000 

1000000 


TABLE IV: Distribution of 1 000 000 random pure states on the different entanglement polytopes. In the first column we list the union of the 
different permuted polytopes of the same type. 























